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t::^- : ABSTRACT 

' We use galaxy groups selected from the Sloan Digital Sky Survey to examine the align- 

ment between the orientation of the central galaxy (defined as the brightest group 
, member) and the distribution of satellite galaxies. By construction, we therefore only 

address the alignment on scales smaller than the halo virial radius. We find a highly 
significant alignment of satellites with the major axis of their central galaxy. This 
. is in qualitative agreement with the recent study of Brainerd (2005), but inconsis- 

JZLh' tent with several previous studies who detected a preferential minor axis alignment. 

The alignment strength in our sample is strongest between red central galaxies and 
^ ■ red satellites. On the contrary, the satellite distribution in systems with a blue central 

^ ' galaxy is consistent with isotropic. We also find that the alignment strength is stronger 

, in more massive haloes and at smaller projected radii from the central galaxy. In ad- 

dition, there is a weak indication that fainter (relative to the central galaxy) satellites 
are more strongly aligned. We present a detailed comparison with previous studies, 
, and discuss the implications of our findings for galaxy formation. 

■ Key words: dark matter - large-scale structure of the universe - galaxies: haloes - 

galaxies: structure - methods: statistical 
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1 INTRODUCTION 

During the hierarchical assembly of dark matter haloes, pro- 
genitor haloes often survive accretion onto a larger system, 
thus giving rise to a population of subhaloes. If the baryonic 
material in the progenitor haloes managed to cool and form 
stars before being accreted, the population of subhaloes will 
give rise to a population of satellite galaxies. Meanwhile, gas 
that cools onto the center of the parent halo gives rise to a 
so-called central galaxy. 

Since satellite galaxies are typically distributed over 
the entire dark matter halo, they are ideally suited as a 
tracer population of the potential well in which they or- 
bit. Consequently, they have been used extensively as dy- 
namical tracers of the dark matter mass distribution sur- 
rounding central galaxies. In addition to providing accurate 
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dynamical masses of the haloes (e.g., Zaritsky et al. 1993, 
1997a; McKay et al. 2002; Brainerd & Specian 2003; van den 
Bosch et al. 2004), the radial trend of the projected velocity 
dispersion of satellite galaxies can also put constraints on 
the radial density distribution of the dark matter (Prada et 
al. 2003). Similar constraints can also come from the mea- 
surement of individual satellite orbits, such as that of the 
Large Magellanic Cloud or the Sagittarius stream in the 
Milky Way halo (e.g., Ibata et al. 2001; Helmi 2004; Kalli- 
vayahl et al. 2006). 

In addition to these kinematics, the spatial distribution 
of satellite galaxies also holds important information. If sub- 
haloes are a fair tracer of the dark matter mass distribution, 
i.e., if they are not spatially biased in any way, the radial and 
angular distribution of satellite galaxies directly reflects the 
projected distribution of the dark matter. If, on the other 
hand, there is a spatial bias, the satellite distribution holds 
important clues regarding the actual assembly history of the 
dark matter haloes. 



2 X. Yang et al. 



Table 1. The observational search for possible alignment between central galaxies and their satellites. 



Attempt 


N systems 


central galaxy type 


Tp/ kpc 


At;|/ km s ^ 




alignment 


LG included 


(1) 


(2) 


(3) 


(4) 


(5) 


(6)' 


(7) 


(8) 


Holmberg (1969) 


218 


nearby spirals 


<50 






minor axis 


YES 


Zaritsky et al. (1997) 


115 


nearby spirals 


<200 


<500 




NO 










[300 500] 


<500 




minor axis 




Sales & Lambas (2004) 


1276 


2dFGRS BIG(*) 


<500 


<160 


<0.16 


minor axis 




Kroupa et al. (2005) 


11 


Milky- Way 


<250 






minor axis 


YES 


Brainerd (2005) 


3292 


SDSS-DR3 BIG 


<700 


<1000 


<0.16 


major axis 






1575 




<500 


<1000 


<0.25 


major axis 






935 




<500 


<1000 


< 0.125 


major axis 




AB (**) (2005b) 


4327 


SDSS-DR4 BIG 


<700 


'<500 


'<0.25 


major axis 




This work (2006) 


24728 


SDSS-DR2 GCG(***) 




<'yvir 




major axis 





Column (1) indicates the attempt ID. Columns (2) (number of central-satellite systems), (3) (the type of central galaxy), (4) (the 
projected centric-distance of the satellite galaxy) (5) (the line-of-sight velocity difference of the satellite-central system), and (6) (the 
luminosity fraction of the satellite-central system) indicate the selection criteria. Column (7) lists the alignment signal obtained from this 
observation. Column (8) lists the status of the Local Group. (*) BIG means Bright Isolated Galaxies, (**) AB: Agustsson & Brainerd 
(2005b) and (***) GCG means Group Central Galaxy. 



Numerical simulations predict that the spatial distribu- 
tion of subhaloes is biased in two distinct ways. First of all, 
the radial distribution of subhaloes is found to be less cen- 
trally concentrated than the dark matter, with a pronounced 
deficit of subhaloes near the center (Ghigna ct al. 1998, 2000; 
Colin et al. 1999; Springel et al. 2001; De Lucia et al. 2004; 
Gao et al. 2004; Diemand, Moore & Stadel 2004; Mao et 
al. 2004; Nagai & Kravtsov 2005; Kang ot al. 2005). A sim- 
ilar prediction has also been obtained with detailed, semi- 
analytical models (Zentner et al. 2005a, but see also Taylor 
& Babul 2004). Somewhat surprisingly, the observed spatial 
distribution of satellite galaxies, especially in clusters, ap- 
pears to bo more centrally concentrated and not as strongly 
anti-biased as this predicted distribution of subhaloes (Carl- 
berg, Yee & EUingson 1997; van den Marel et al. 2000; Lin, 
Mohr & Stanford 2004; van den Bosch et al. 2005; Yang et 
al. 2005c; Chen et al. 2005). A possible explanation, which 
needs to be explored in detail, is that the addition of baryons 
makes the subhaloes more resilient to tidal disruption. 

In addition to this radial anti-bias, numerical simula- 
tions have also suggested an angular bias. In particular, 
numerous studies have shown that dark matter haloes in 
dissipationless simulations have anisotropic distributions of 
subhaloes that are aligned with their major axis (Knebe et 
al. 2004; Zentner et al. 2005b; Wang et al. 2005; Libeskind 
et al. 2005). This anisotropy mainly owes to a prefered di- 
rection of satellite accretion along large scale filaments (Tor- 
men 1997; Vitvitska ct al. 2002; Knebe et al. 2004; Aubert, 
Pichon & Colombi 2004: Zoutner et al. 2005b; Wang et 
al. 2005). Since the orientation of the halo itself is largely 
governed by the directionality of its mass accretion (e.g., 
van Haarlem & van dc Weygacrt 1993; Tormen 1997), this 
naturally explains the alignment of the subhaloes with the 
major axis of the parent halo. It is important to distinguish 
here between simple angular anisotropy of the subhalo dis- 
tribution, and a true angular bias. Even in the absence of 
any spatial bias, any non-sphericity of dark matter haloes 
will result in a projected, angular anisotropy of the subhalo 
distribution, unless the halo is seen along its symmetry axis. 
Both Wang et al. (2005) and Agustsson & Brainerd (2005a) 



have shown that the halo's non-sphericity is the main cause 
of the theoretical angular anisotropy, but that there is a 
weak indication for some additional angular bias with a pre- 
ferred alignment along the halo's major axis. 

If the orientations of central galaxies are somehow 
aligned with their dark matter haloes, this anisotropy should 
result in an observable correlation between the distribu- 
tion of satellite galaxies and the orientation of their cen- 
tral galax;y. Numerical simulations suggest that the angular 
momenta of dark matter haloes are typically aligned with 
their minor axes (e.g.. Warren ct al. 1992; Dubinski 1992; 
Porciani, Dekel & Hoffman 2002a; Bailin & Steinmetz 2005; 
Faltenbacher et al. 2005; Wang et al. 2005). If the angular 
momentum vector of the baryonic material is well aligned 
with that of the dark matter, one would thus naively expect 
the spin axes of disk galaxies to be aligned with the minor 
axes of their host haloes, and the satellite galaxies to be 
preferentially oriented along the disk's major axis. However, 
detailed hydro-dynamical simulations reveal a more compli- 
cated picture. First of all, even in the absence of cooling 
the spin axes of the baryons and the dark matter are only 
poorly aligned, with a median misalignment angle of about 
~ 20° to ~ 30° (van den Bosch et al. 2002; Chen, Jing & 
Yoshikaw 2003; Sharma & Steinmetz 2005). Furthermore, 
in simulations of disk galaxy formation that include cooling, 
the orientation of the disk spin axis is found to be virtu- 
ally uncorrelated with the original (i.e., in the absence of 
baryons) minor axis of the halo (Bailin et al. 2005). The 
formation of the disk modifies the shape and orientation 
of the inner halo, but leaves the outer halo largely intact 
(Kazaiitzidis et al. 2004; Bailin et al. 2005). Consequently, 
the disk spin axis is well aligned with the halo minor axis 
in the inner halo (r <0.1rvir), but is basically uncorrelated 
with the minor axis at larger halo-ccntric radii. If correct, 
this would predict basically no alignment between the orien- 
tation of the central disk galaxy and the distribution of its 
satellites (most of which lie at relatively large halo-centric 
radii) . 

The observational search for a possible alignment be- 
tween central galaxies and their satellites has a long and 



The Alignment between Satellites and their Central Galaxy 3 



confusing history. Holmberg (1969) studied the distribu- 
tion of sateUite galaxies around isolated disk galaxies, and 
found them to lie preferentially along the minor axis of 
disk galaxies. Holmberg's study was restricted to projected 
satellite-central distances of Vp < 50kpc. Subsequent studies, 
however, were unable to confirm this so-called "Holmberg- 
effect" (Hawley & Peebles 1975; Sharp, Lin & White 1979; 
MacGillivray et al. 1982). Zaritsky et al. (1997b) were also 
unable to detect any significant alignment for Vp < 200 kpc, 
but they did detect a preferred minor axis alignment for 
separations in the range 300 kpc<rp <500 kpc. Note that 
this implies an alignment on scales larger than the typical 
virial radii of the haloes hosting these isolated disk galax- 
ies. This large-scale (rp < 500 kpc) alignment has recently 
been confirmed by Sales & Lambas (2004), using data from 
the two-Degree Field Galaxy Redshift Survey, but only for 
host-satellite pairs with a line-of-sight velocity difference of 
|A«| < leOkms"^ Our own Milky- Way also reveals a Holm- 
berg effect, in that the 11 innermost MW satellites (with 
MW distances < 250 kpc) show a pronounced planar dis- 
tribution oriented close to perpendicular to the MW disk 
(Lynden-Bell 1982; Majewski 1994; Kroupa, Theis & Boily 
2005). Completely opposite to all theses results, Brainord 
(2005) and Agustsson & Brainerd (2005) recently found 
that, in the Sloan Digital Sky Survey (SDSS), the distribu- 
tion of satellite galaxies with Vp < 100 kpc is strongly aligned 
with the major axis of the disk host galaxy. As a summary, 
we list in Table 1. the main attempts in searching for the 
alignment signal between the central galaxies and their satel- 
lites. 

Clearly, this lack of agreement calls for a more in-depth 
study. In this paper we investigate the alignment between 
satellite galaxies and their host galaxies using data from the 
SDSS. Our approach, however, differs substantially from all 
previous studies. First of all, previous studies only focused 
on relatively isolated disk galaxies, which has the disadvan- 
tages that it drastically reduces the sample size, and that 
one only selects haloes in relatively low-density environ- 
ments. In addition, satellite galaxies were always selected 
in a fixed metric aperture centered on the host galaxy. For 
low luminosity hosts, which reside in low mass haloes, this 
metric is often much larger than the expected virial radius 
of the host halo. In this paper we study the host-satellite 
alignment using a large sample of galaxy groups. No iso- 
lation criteria are applied, which allows us to (i) achieve 
much better statistics, and (ii) to investigate how the align- 
ment strength depends on various properties of the host 
halo, the host galaxy, and the satellite galaxies. In addi- 
tion, wo only focus on satellites that are located within the 
host halo's virial radius with projected satellite-central dis- 
tances Tp < Tvir and satellite-central line-of-sight velocity 
differences |Ai;| < Vvir, where rvir and Wvir are the virial ra- 
dius and virial velocity dispersion of the host dark matter 
halo, respectively. In other words, we select satellites using 
a variable aperture size that is motivated by the mass of the 
host halo (i.e., the galaxy group). This has the important ad- 
vantage that we can clearly separate small-scale alignment 
(r < Tvir) from large-scale alignment (r > rvir). Depending 
on how the shapes and angular momenta of dark matter 
haloes are oriented with respect to their surrounding large 
scale structure, the alignment on these two different scales 
may well be very different (see e.g., Barnes & Efstathiou 



1987; Porciani et al. 2002a,b; Navarro, Abadi & Steinmetz 
2004; Bailin & Steinmetz 2005; Trujillo, Carretero & Patiri 
2005). 

This paper is organized as follows. In Section 2 we 
present the data and our methodology. Section 3 presents a 
careful analysis of the alignment between the orientation of 
central galaxies and the distribution of their satellite galax- 
ies. In particular, we show how the alignment strength de- 
pends on the luminosity and color of the galaxies, and on the 
mass of the dark matter haloes. In Section 4 we present a de- 
tailed comparison with previous studies. Finally, in Section 5 
we summarize our results, and discuss their implications. 



2 DATA AND ANALYSIS 
2.1 Galaxy groups 

In order to address the possible alignment between satel- 
lite galaxies and the central galaxy of their dark matter 
parent halo we use the SDSS galaxy group catalogue of 
Weinmann et al. (2005)^. This catalogue was constructed 
from the New York University Value-Added Galaxy Cata- 
logue (NYU-VAGC; Blanton et al. 2005) ^ using the halo- 
based group finder developed by Yang et al. (2005a, here- 
after YMBJ). The NYU-VAGC is based on the SDSS Data 
Release 2 (Abazajian et al. 2004), but with an independent 
set of significantly improved reductions. We only consider 
the galaxies with redshifts in the range 0.01 < z < 0.2 and 
with a redshift completeness c > 0.7 ^, resulting in a sample 
of 184,425 galaxies with a sky coverage of ~ 1950 deg^. 

In brief, the YMBJ group finder works as follows. First 
potential group centers are identified using a Friends-Of- 
Friends (FOF) algorithm or an isolation criterion. Next, the 
total group luminosity is estimated which is converted into 
an estimate for the group mass using an assumed mass-to- 
light ratio. From this mass estimate, the radius and velocity 
dispersion of the corresponding dark matter halo are esti- 
mated using the virial equations, which in turn are used 
to select group members in redshift space. This method is 
iterated until group memberships converge. Detailed tests 
with mock galaxy redshift surveys have shown that this 
group finder recovers groups with an average completeness 
of ^ 90% and with an interloper fraction that is smaller than 
~ 20%. The resulting group catalogue is insensitive to the 
initial assumption regarding the mass-to-light ratios, and 
the group finder is more successful than the conventional 
FOF method in associating galaxies according to their com- 
mon dark matter haloes (see YMBJ for details). 

Following Yang et al. (2005b) we use the group luminos- 
ity to assign masses to our groups. The motivation behind 
this is that one naturally expects the group luminosity to 
be strongly correlated with halo mass (albeit with a certain 
amount of scatter). For each group we determine the number 

^ In this paper, we refer the brightest member in each group as 
the central galaxy, while all other members as satellite galaxies. 
^ http:/ /wassup. physics. nyu.edu/vagc/#:download 

Because of the survey selection effects (e.g., fiber collisions, 
etc.), not all galaxies in the photometric catalogue are spec- 
troscopically observed and thus their spectroscopic redshifts are 
measured. 
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Figure 1. The normalized probability distribution of the angle 
6 between the orientation of the major axis of the central group 
galaxy and the direction of each satellite as measured from the 
central galaxy. These results have been obtained from the SDSS 
group catalogue discussed in the text, were we have excluded 
those groups for which the projected ellipticity of the central 
galaxy is less than 0.2. This leaves a grand total of 24, 728 unique 
central-satellite pairs. The normalization and crrorbars are com- 
puted from 100 random samples in which we have randomized 
the orientation of all central galaxies (see text for details). Note 
that /pairs > 1 for 6 < 35° indicating that the satellite galaxies 
are preferentially distributed along the major axis of their central 
galaxy. This is also evident from the average value of 6, and its 
error, which are indicated in the upper-right corner. Note that an 
isotropic satellite distribution corresponds to (0) = 45° . 



density of all groups brighter than the group in considera- 
tion, using a common, empirically calibrated definition of 
group luminosity. FYom the halo mass function correspond- 
ing to a ACDM concordance cosmology with flm = 0.3, 
fiA = 0.7, h = Ho/{WO kms"^ Mpc"^) = 0.7 and as = 0.9 
we then find the mass for which the more massive haloes 
have the same number density. Although this has the down- 
side that it depends on cosmology, as shown in Weinmann 
et al. (2005), this method yields masses that are more accu- 
rate than those based on the more traditional line-of-sight 
velocity dispersion of the group members (see Weinmann et 
al. 2005). In addition, it is straightforward to convert the 
masses derived here to any other cosmology (see Yang et 
al. 2005b). 

Applying our group finder to the sample of 184, 425 
galaxies in the NYU-VAGC described above yields a group 
catalogue of 53, 229 systems with an estimated mass. These 
groups contain a total of 92, 315 galaxies. The majority of 
the groups (37, 216 systems) contain only a single member, 
while there are 9220 binary systems, 3073 triplet systems. 



and 3720 systems with four members or morc^. In what fol- 
lows, we use this group catalogue to examine the alignment 
between the orientation of the central galaxy, defined as the 
brightest group member, and the distribution of satellite 
galaxies. Note that we have a total of 39, 086 unique central- 
satellite pairs, which is an order of magnitude larger than in 
any previous study. 



2.2 Methodology 

In order to quantify the distribution of satellite galaxies in 
groups relative to the orientations of their central galax- 
ies we follow Brainerd (2005) and compute the distribution 
function, P{0), where is the angle on the sky between 
the major axis of the 25 magnitudes per square arcsecond 
isophote in the r-band of the central group galaxy and the 
direction of a satellite relative to the central galaxy. We re- 
strict to the range [0°, 90°], where = 0° (90°) implies that 
the satellite lies along the major (minor) axis of the central 
galaxy. The orientation of the central galaxy is based on the 
isophotal position angle in the r-band, as given in the SDSS 
Data Release 2 (Abazajian et al. 2004). 

As mentioned above, the analysis of Brainerd (2005) 
focused on relatively isolated systems with late-type central 
galaxies. Our analysis is different in that we consider galaxy 
groups of various properties. In practice, we start with a 
group sample and count the total number of central-satellite 
pairs, N{9), for a number of bins in 9. Next we construct 
100 random samples in which we randomize the orientation 
of all central galaxies, and we compute {Nii{9)), the average 
number of central-satellite pairs as function of 9. Note that 
this ensures that the random samples have exactly the same 
selection effects as the real sample, so that any significant 
difference between A''(^) and Nr{9) reflects a genuine align- 
ment between the orientation of the central galaxies and the 
distribution of satellite galaxies. 

To quantify the strength of any possible alignment we 
define the normalized pair count. 



/pairs(e) = N{9)/(Nr{9)) , 



(1) 



Note that /pairs (S) = 1 in the absence of any alignment. We 
use a r{9) / {N ii{9)) , where aR{d) is the standard deviation 
of Nii{9) obtained from the 100 random samples, to assess 
the significance of the deviation of /pairs (^) from unity. In 
addition to this normalized pair count, wo also compute the 
average angle (9) . In the absence of any alignment (9) = 45° , 
however, (9) = 45° does not mean an isotropic distribution. 
The significance of any alignment can be expressed in terms 
of {9} and ae, the variance in {9)r as obtained from the 100 
random samples. 

Finally, since the accuracy with which 9 can be mea- 
sured scales with the projected ellipticity of the central 
galaxy, we only consider groups for which the ellipticity of 
the central galaxy, e > 0.2. Here e is defined as one mi- 
nus the ratio between the minor and major axes of the 25 
magnitudes per square arcsecond isophote in the r-band of 
the image of the central galaxy. This ellipticity constraints 



The group catalogue is publicly available 
http://www.astro.umass.edu/~xhyaiig/Group.html 
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Figure 2. The same as Fig. 1 but for different subsarnples of central and satellite galaxies. In the upper panels, we show /pairs 
for groups with a different cUipticity, e, of the central galaxy, as indicated. Note that groups with a strongly elongated central galaxy 
(0.6 < e < 0.8) are consistent with a perfectly isotropic distribution of satellites. As we argue in the text, and show in Fig 3, this owes to 
the fact that strongly elongated systems are mainly blue, late type disk galaxies, which show no significant alignment. The lower panels 
show how /pairs (^) depends on the luminosities of the satellite galaxies, L^, expressed in units of the luminosity of their central galaxy, 
Lc- There is a clear indication that fainter satellites are more strongly aligned. 



brings the total number of unique central-galaxy pairs to 
24, 728 



3 RESULTS 

Fig. 1 shows /pairs (0) foi all groups in our SDSS group cat- 
alogue with inferred halo masses of M > 10^^ ft-^Mo and 
with central galaxies that have e > 0.2. Note the pronounced 
enhancement of pairs with small 6, implying that satellite 
galaxies are preferentially distributed along the major axes 
of their central galax;ies. This is also evident from the fact 
that {&) = 42.2° ± 0.2°, which deviates from the case of no 
alignment (i.e., (6') = 45.0°) by 14ct! 

Since the accuracy of the orientation angle of a central 
galaxy is smaller for central galaxies that appear rounder, 
the strength of the alignment may be diluted due to cen- 
tral galaxies with a small ellipticity, e. In order to address 
the impact of e on the strength of the alignment signal, the 
upper panels of Fig. 2 show /pairs (^) for groups with cen- 
tral galaxies with different oUipticities, as indicated. Note 
that the alignment strength is weakest for the sample with 
the highest ellipticities (0.6 < e < 0.8). This is surprising 
since one would expect the orientation angle of these cen- 
tral galaxies to be the most accurate. However, as we will see 
in Section 3.1, the strength of the alignment is significantly 
weaker for systems with blue, late-type central galaxies than 



for systems with red, early-type central galaxies. The depen- 
dence on e found here is simply due to the fact that central 
galaxies with e > 0.6 are dominated by blue, late-type, disk 
galaxies. In what follows we always use all groups with cen- 
tral galaxies with e > 0.2. 

3.1 Dependence on Galaxy Properties 

Since our group catalogue contains a large number (~ 
39, 000) of central-satellite pairs, it allows us to study how 
the alignment depends on various properties of the central 
and satellite galaxies. We start by examining the depen- 
dence on the luminosities of the satellite galaxies. The lower 
panels of Fig 2 show /pairs (^) for a number subsamples of 
satellite galaxies that are selected based on their luminosi- 
ties, Ls, relative to the lumiuositics of their central galaxies, 
Lc. There is a weak trend that fainter satellite galaxies are 
more strongly aligned with the orientation of their central 
galaxy than brighter satellite galaxies. 

Next, we consider the dependence on the color of the 
satellite galaxies. We separate galaxies into two subsamples 
according to their '^'^{g — r) colors, which corresponds to 
the [g — r) color fc-corrected to redshift z = 0.1. We call 
galaxies with "'^{g — r) < 0.83 'blue' and galaxies with 
"'^{g—r) > 0.83 'red'. This value of 0.83 roughly corresponds 
to the bimodality scale in the color-magnitude relation (see 
Weinmann et al. 2005). 
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Figure 3. The same as Fig. 1, but for different subsamples of 
liosts and satellites, selected according to their "'^(g — r) color. 
See text for discussion. 



Figure 4. Same as Fig 3, except that here we split the sample ac- 
cording to the colors of both the central and the satellite galaxies, 
as indicated. 



The upper panels of Fig. 3 show /pairs {8) for blue and 
red satellite galaxies, while the lower panels show /pairs (5) 
for groups with blue and red central galaxies. Note that 
there is a remarkably strong dependence on the color of both 
the central galaxies and the satellite galaxies. In particular, 
satellite galaxies in groups with a blue, central galaxy are 
consistent with a perfectly isotropic distribution; there is no 
sign of any significant alignment {{6) = 44.5° ±0.5°). On the 
contrary, groups with a red central galaxy show a very pro- 
nounced, major-axis alignment with (6) — 41.5° ±0.2°. In 
addition, red satellites show a significantly stronger major 
axis alignment than blue satellites. 

As shown in Weinmann et al. (2005), haloes with a cen- 
tral red galaxy have a significantly larger fraction of red 
satellites than a halo of the same mass, but with a blue cen- 
tral galaxy. This so-called 'galactic conformity' implies that 
the upper and lower panels are not independent. In Fig. 4 
we therefore examine how /pairs (S) depends on the colors of 
both the central galaxy and the satellites. As can be seen, 
systems with a blue central galaxy show no significant align- 
ment, neither with their blue satellites nor with their red 
satellites. Systems with a red central galaxy, however, show 
a very pronounced alignment, which is significantly stronger 
for red satellites than it is for blue satellites. Since redder col- 
ors typically indicate older stellar populations, these results 
suggest that a significant alignment between the orienta- 
tion of central galaxies and the distribution of their satellite 
galaxies only exists in haloes with a relatively old stellar 
population. Clearly, such a correlation between the align- 
ment strength and the age of the stellar population must 
hold some interesting clues regarding galaxy formation. We 
return to this in Section 5. 



3.2 Halo Mass Dependence 

It is interesting to examine whether the alignment strength 
also depends on halo mass. Since our group catalogue covers 
a large range in halo masses, we can address this question in 
some detail. The upper panels in Fig 5 show the results for 
groups in 3 mass bins, as indicated. There is a clear mass- 
dependence, in the sense that the alignment is stronger for 
more massive groups. 

Since more massive haloes contain a larger fraction of 
red galaxies (e.g., Weinmann et al. 2005), and given that red 
galaxies show a much more pronounced alignment than blue 
galaxies, this mass dependence may simply reflect the color 
dependence shown in the previous section. To test this, the 
panels in the middle and lower row of Fig 5 show /pairs (^) for 
groups in the same three mass bins, but considering groups 
with blue or red central galaxies separately. This shows that 
haloes with a blue central galaxies show no significant align- 
ment, independent of their mass. Haloes with a red cen- 
tral galaxy, however, do show a significant alignment, with 
a strength that increases with increasing halo mass. Thus, 
there is a genuine mass dependence, but only for haloes with 
red central galaxies. 

3.3 Radial Dependence 

Finally, we examine whether the alignment depends on the 
group-centric distance of satellite galaxies. For each central- 
satellite pair we compute the projected separation, r, in 
units of the virial radius, rvir, of the corresponding halo. 
Fig. 6 plots {6) as function of r/rvir- An isotropic satellite 
distribution will have (6) — 45° , indicated by the horizontal 
line, while values smaller (larger) than 45° indicate a pre- 
ferred alignment with the major (minor) axis of the central 
galaxy. As before, the errorbars are obtained from the 100 
randomizations of the orientations of the major axes of the 
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Figure 5. The same as Fig 1, but for satellite-central galaxy pairs in haloes of different mass bins as indicated at the top of the panels. 
In the upper panels we show the results using all groups with e > 0.2, independent of the color of the central galaxies. Note that there 
is a weak trend of increasing alignment strength with increasing halo mass. Panels in the middle and lower rows correspond to haloes 
in the same mass ranges but with only blue and red central galaxies, respectively. Note that haloes with blue, central galaxies show no 
significant alignment of their satellite distribution with the orientation of the central galaxy, independent of halo mass. Haloes with a 
red, central galaxy on the other hand, always reveal a major axis alignment, with a strength that increases with halo mass. 



central galaxies. The upper, left-hand panel shows the re- 
sults for all groups with central galaxies with e > 0.2. There 
is a clear radial trend, in that satellites at smaller, projected 
distances from their central galaxy are more strongly aligned 
with its major axis. If we split the sample according to the 
satellite luminosity, Ls , relative to that of the central galaxy, 
Lc, there is a pronounced difference: fainter satellites do not 
show a significant radial dependence, while satellites with 
Ls > 0.3I/C show a very strong radial trend. 

The lower panels of Fig. 6 show (6) as function of r /rvir 
for three different bins in halo mass, as indicated. Haloes 
with M < lO^^/i~^M0 reveal a significant trend of decreasing 
{9} (i.e., a stronger major-axis alignment) with decreasing 
radius. In more massive haloes, however, there is no signif- 



icant radial trend. Instead, in these haloes the major axis 
alignment is extremely strong at all projected radii. 



Fig. 7 shows how {9){r/r^ii) depends on the colors of 
the satellites and their central galaxies. Systems with a blue 
central galaxy only show a weak (~ 3a") major axis align- 
ment with satellites (both red and blue) at r < 0.2r-vir. 
The distribution of satellites at larger projected radii is per- 
fectly consistent with isotropic. Systems with a red, central 
galaxy reveal a weak, but significant, radial trend of decreas- 
ing alignment strength with increasing radius. This is most 
pronounced for the red satellites, while the blue satellites 
are consistent (within the errors) with a constant alignment 
strength at all projected radii. 
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Figure 6. The average angle 6 as function of tlie projected radius between the sateUite galaxy and the central group galaxy, r, expressed 
in units of the group's virial radius, r^ir- The upper left-hand panel shows the result for all groups in which the central galaxy has an 
elliptieity e > 0.2. The upper panels in the middle and to the right correspond to central-satellite pairs for Ls/Lc falls in the range 
indicated. The lower panels show the results for three subsamplcs selected according to the mass of the groups, again as indicated. The 
thin, horizontal line indicates (9) = 45°, which corresponds to an isotropic distribution. The errorbars arc computed from 100 random 
samples in which we randomized the orientation of all central galaxies. Overall, the major axis alignment strength is stronger at smaller 
projected radii. 



4 COMPARISON WITH PREVIOUS STUDIES 

As shown above, wc detect a significant alignment of satel- 
lite gala^xies with the major axis of their central host galaxy. 
This is in qualitative agreement with the recent studies of 
Brainerd (2005, hereafter BOS) and Agustsson & Brainerd 
(2005b), but in strong disagreement with Holmberg (1969), 
Zaritsky ot al. (1997b, hereafter ZSFW97) and Sales & Lam- 
bas (2004, hereafter SL04). 

First of all, given that many studies have been unable to 
reproduce the results of Holmberg (1969), and given that he 
used a sample consisting of only 58 hosts and 218 satellites, 
we argue that Holmberg's results are probably an unfortu- 
nate outcome of the small sample size. Secondly, our results 
are not necessarily inconsistent with those of ZSFW97, who 
only detected a significant minor axis alignment at relatively 
large projected radii (300 kpc <r< 500 kpc) . This is larger 
than the typical virial radius, Tvir, expected for the isolated 
disk galaxies used in their study. Since we have only fo- 
cused on the alignment at scales r < r^ir, our results and 
theirs are not mutually exclusive. For r < 200 kpc ZSFW97 
did not find any indication for a significant satellite align- 
ment. Since they only focused on isolated disk galaxies, this 
is in agreement with the isotropic distribution of satellites 
in systems with a blue central galaxy presented here. Note 
that ZSFW97 only had a sample consisting of 115 satellites 
(around 69 host galaxies). 



The results of SL04, based on the 2dFGRS, are more dif- 
ficult to explain in light of our findings. In particular, SL04 
also investigated how their alignment strength correlates 
with the properties of the central galaxies. In agreement 
with our results, they find that the satellite distributions 
around blue centrals arc consistent with being isotropic, 
while satellites around red centrals show a strong alignment 
effect. However, contrary to the major axis alignment found 
here, SL04 detected a minor axis alignment. The fact that 
the same trends are detected, but in the opposite direction, 
is suggestive of an error in the computation of 9, and we 
have performed a number of tests to investigate this possi- 
bility. Unfortunately, the major axis position angles of the 
2dFGRS galaxies are defined as "measured in degrees clock- 
wise from East to West"^. This description is ambiguous as 
it is unclear whether an angle of 45° corresponds to North- 
east (as would be the case if the astronomical convention is 
used) or to Southeast. To test this, we cross-correlate the 
2dFGRS with our SDSS sample, and compare the orienta- 
tion angles provided by both catalogues. This comparison 
indicates that the 2dFGRS orientation angles are measured 
from East through South. As it turns out, SL04 interpreted 
the orientation angles as running from East through North 
(Laura Sales and Diego Lambas, private communication). 



^ http: / /www.mso.anu.edu.au/2dFGRS/ 
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Consequently, what they call a minor axis alignment is in 
fact a major axis alignment. In retrospect, the SL04 results 
are thus in qualitative agreement with BOS and with the 
results presented here. This has been confirmed by tests 
performed by the authors (Laura Sales and Diego Lambas, 
private communication). 

Finally, it is worth pointing out that our qualitative 
agreement with the study of B05, who also used SDSS data, 
is not entirely trivial. While B05 focused on relatively iso- 
lated central galaxies, we consider all galaxy systems, from 
poor groups to rich clusters, selected with our group finder. 
Furthermore, BOS only select systems in which all satellites 
are significantly fainter than the central galaxy (similar to 
ZSFW97 and SL04). In fact, the selection criteria are so re- 
strictive, that although BOS starts out with a larger SDSS 
sample than used here (SDSS Data Release 3 versus SDSS 
Data Release 2), she ends up with less than 3300 satellites, 
almost an order of magnitude fewer than in our case. Given 
these dramatic differences, it is therefore not obvious that 
the results have to be compatible. We have verified, however, 
that if we select, from our group catalogue, those groups that 
obey the selection criteria of BOS, we obtain results that are 
of slightly stronger alignment signal than those of BOS. This 
discrepancy may partly due to the fact that the color distri- 
bution of the remaining central galaxies are very similar to 
the original ones, unlike those in BOS which are mostly spiral 
galaxies and, in our investigations, have smaller alignment 
signal. 

To summarize, BOS and SL04 both have obtained results 
that are in qualitative agreement with the results presented 
here: Satellites around red hosts are aligned with the host's 
major axis, while satellites around blue hosts have an angu- 
lar distribution that is consistent with isotropic. As we have 
argued, these results are not necessarily inconsistent with 
those of ZSFW97. The only study that is in clear disagree- 
ment with our results is that of Holmberg (1969). However, 
given his relatively small sample size, this discrepancy is 
not very significant. Nevertheless, one important exception 
to the general picture obtained here exists, namely our own 
MW. As discussed in Section 1, the MW satellites reveal a 
pronounced planar distribution that is oriented close to per- 
pendicular to the disk. However, as we argue below, this is 
not necessarily inconsistent with our results. 



5 DISCUSSION &: CONCLUSIONS 

Using galaxy groups selected from the SDSS, we have ex- 
amined the alignment between the orientation of the cen- 
tral galaxy (defined as the brightest group member) and 
the distribution of its satellites. Overall we find an excess 
of satellites along the major axis, and a deficiency along the 
minor axis, compared to an isotropic distribution. The align- 
ment strength in our sample is strongest between red central 
galaxies and red satellites. On the contrary, the satellite dis- 
tribution in systems with a blue central galaxy is perfectly 
consistent with isotropic. We also find that the alignment 
strength is stronger in more massive haloes and at smaller 
projected radii from the central galaxy. In addition, there is 
a weak indication that fainter (relative to the central galaxy) 
satellites are more strongly aligned. 

Two conditions must be satisfied in order to produce 
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Figure 7. Same as Fig. 6 except that we show the results for 
different subsamples selected according to the "'^(g — r) colors 
of the central and satellite galaxies, as indicated. See text for a 
detailed discussion. 



the alignment observed here. First of all, the distribution of 
satellite galaxies in groups must be aspherical, and secondly, 
the orientation of central galaxies must be aligned with the 
distribution of satellite galaxies. Cosmological A'^-body simu- 
lations of CDM models have demonstrated clearly that CDM 
haloes are not spherical. The typical minor-to-major axis ra- 
tio is ~ 0.6, with a relatively large dispersion (e.g. Bullock 
2002; Jing & Suto 2002). This is also supported by recent 
weak-lensing data (Hoekstra, Yee & Gladders 2004). As we 
argued in Section 1, simulations suggest that the angular 
distribution of dark matter subhaloes, which are expected 
to host satellite galaxies, is in reasonable agreement with 
that of the dark matter. Therefore, it seems reasonable to 
expect that the first condition is fulfilled. Explaining the 
anisotropy as reflecting the non-sphericity of the dark mat- 
ter haloes is also consistent with our flnding that the align- 
ment strength increases with halo mass; after all, numerical 
simulations have shown that more massive haloes are less 
spherical (Warren et al. 1992; Bullock 2002; Jing & Suto 
2002; Bailin & Steinmetz 2005; Kasun & Evrard 200S). 

In order for the second condition to be fulfilled as well, 
the central galaxy must be somehow aligned with the prin- 
cipal axes of the mass distribution of its host halo. Here it is 
important to distinguish between disk galaxies, whose orien- 
tation is governed by their angular momentum vector, and 
spheroidal galaxies, whose orientation is somehow related 
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to its formation history (typically thought to be merger- 
driven). 

In the standard picture of disk formation (e.g., Fall & 
Efstathiou 1980; Mo, Mao & White 1998), one assumes that 
baryons and dark matter have identical distributions of spe- 
cific angular momentum (due to tidal torques from the cos- 
mological density field) , and that the baryons conserve their 
specific angular momentum when cooling to form a centrifu- 
gally supported disk. Since simulations have shown that the 
spin axis of dark matter haloes is well aligned with the halo's 
minor axis, this simple picture predicts that the disk spin 
axis should be parallel to the minor axis of the halo, and thus 
that the distribution of satellites is aligned with the disk 
major axis. Somewhat surprisingly, disk galaxies (which are 
typically blue) are exactly the subsample that do not seem 
to reveal a significant alignment with their satellites. 

Detailed hydrodynamical simulations, however, have 
shown that the spin axes of the baryons and the dark mat- 
ter (in the absence of cooling) are only poorly aligned (van 
den Bosch et al. 2002; Chen et al. 2003; Sharma & Stein- 
metz 2005). In addition, if cooling is included in the simu- 
lations, the resulting disks are found to have spin axes that 
are very poorly aligned with the original (i.e., in the absence 
of baryons) minor axis of the halo (Bailin et al. 2005). This 
suggests that the (direction) of the angular momentum of 
the baryons is not well conserved during the disk formation 
process, and thus predicts that there is little if any alignment 
between the orientation of the disk and that of its satellite 
distribution. However, there may still be a clear alignment 
between the distribution of the satellites and the principal 
axes of the dark matter halo. This picture not only explains 
the lack of a significant alignment when stacking many disk 
galaxies, but also the pronounced alignment found for the 
MW system: it only requires that the major axis of the MW 
halo happens to be oriented along the spin axis of the MW 
disk. However, in a recent merger-driven disk formation the- 
ory of Robertson et al. (2005) that the disk galaxies can be 
produced through high angular momentum accretion of gas 
rich progenitors, the satellite galaxies will be preferentially 
aligned with the disk plane. Note also that the MW satellites 
galaxies have an average luminosity ratio that is significantly 
smaller than our central-satellite systems. Given this, for the 
MW case, as discussed in Kroupa, Theis & Boily (2005), 
another plausible explanation of the MW (dwarf) satellites' 
distribution is that, if most of the dwarves are not of dark 
matter dominated, but stem from one initial gas-rich parent 
satellite on an eccentric near polar orbit that interacted with 
the young MW, forming tidal arms semi-periodically as its 
orbit shrank, this signal shall be naturally expected. 

Perhaps somewhat surprisingly, we find the strongest 
alignment between red central galaxies and red satellites. 
If the orientation of a red (early-type) central galaxy is a 
reflection of the orbital angular momentum of the progen- 
itors that merged during its formation, one might naively 
expect that this orientation is also aligned with the major 
axis of the halo. After all, the orientation of the halo itself 
is largely governed by the directionality of its mass accre- 
tion (e.g., van Haarlem & van de Weygaert 1993; Tormen 
1997). It is less clear, however, why the alignment should be 
so much stronger for red satellites than for blue satellites. 
In the standard picture of galaxy formation, once a satellite 
galaxy has been accreted by a bigger system, its gas reser- 



voir is stripped, resulting in a fairly quick truncation of star 
formation; the galaxy will become red. In this simple pic- 
ture, one thus expects the color of a satellite galaxy to be a 
reflection of the time since it was accreted. If the orientation 
of a halo, and its population of galaxies, changes as function 
of time with respect to the large scale matter distribution, 
one could envision that those satellites that were accreted at 
around the same time when the central galaxy formed, show 
a more pronounced alignment than those satellites that have 
only recently been accreted. As shown by Bailin & Steinmetz 
(2004), most haloes reveal some slow figure rotation, with 
an amplitude that can cause a directional change of more 
than 90° within a Hubble time. Figure rotation is thus a po- 
tential explanation for the satellite-color dependence of the 
alignment strength. 

Cosmological A'^-body simulations also show that more 
meissive dark matter haloes in general have more elongated 
structures, and that the iso-dcnsity contours of the dark 
matter distributions are strongly aligned in the inner part 
of a halo (e.g. Jing et al. 1995; Jing & Suto 2002; Bailin & 
Steinmetz 2005). Thus, if the central galaxy in such a halo 
aligns with the inner part of the halo, and if the distribution 
of satellite galaxies traces the dark matter distribution, the 
major axis of the central galgixy is expected to align with 
the satellite distribution, and the alignment is expected to 
be stronger on smaller radii. This is in agreement with what 
we found in this paper. 

In a recent study, Augustsson & Brainerd (2005a) used 
the GIF simulations (e.g. Kauffmann 1999) to interpret the 
observed alignment signal between the host and satellite 
galaxies in the context of structure formation. They found 
that the alignment between the satellites and the host halo 
major axis is much stronger than that between the host and 
satellite galaxies in the observations. This would imply that 
the host galaxies are on some degree mis-aligned with the 
host halos. Meanwhile a related interesting result was re- 
cently obtained by Mandelbaum et al. (2006a;b), who used 
the galaxy-galaxy weak lensing signals in the SDSS observa- 
tions to constrain the ellipticity of dark matter haloes. They 
found that, for spirals (lens), the ellipticity of halo and light 
is anti-aligned on a l-2a level, while for ellipticals (lens), the 
ellipticity of halo relative to light, fh = ehaio/eught, increase 
with luminosity. Apparently, these findings axe helpful to in- 
terpret the overall central-satellite galaxy alignment signals 
we obtained in this paper, especially the halo mass, lumin- 
sity and type dependences. 

The discussion presented above shows that the results 
obtained in the present paper are in qualitative agreement 
with naive, theoretical expectations. In order to compare 
our results to theory in a more quantitative manner, one 
has to understand in more detail how central galaxies in 
different halos form, how the formation processes aifect the 
orientations of the central galaxies relative to the haloes, and 
how satellite galaxies trace the mass distribution within dark 
matter halos. All these issues require detailed rmmerical sim- 
ulations of galaxy formation in the cosmic density field. In 
addition, to study the color dependence of the alignment, 
some basic treatment of star formation will be necessary. 
Semi-analytical techniques combined with high-resolution 
numerical simulations may be particularly suited for this 
purpose. We wiU return to these issues in a forthcoming pa- 
per (Kang et al., in preparation). 
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